Energy Budget Research Articles

Seasonality of spectral radiative fluxes and optical properties of Arctic sea ice during the spring–summer transition

The reflection, absorption, and transmittance of shortwave solar radiation by sea ice play crucial roles in physical and biological processes in the ice-covered Arctic Ocean and atmosphere. These sea-ice optical properties, particularly during the melt season, significantly impact energy fluxes within and the total energy budget of the coupled atmosphere-ice-ocean system. We analyzed data from autonomous drifting stations to investigate the seasonal evolution of the spectral albedo, transmittance, and absorptivity for different sea-ice, snow, and surface conditions measured during the MOSAiC expedition in 2019–2020. The spatial variability of these properties was small during spring and increased strongly after melt onset on May 26, 2020, when liquid water content on the surface increased, largely accounting for the enhanced variability. The temporal evolution of surface albedo and sea-ice transmittance was mostly event-driven, thus containing episodic elements. Melt ponds reduced the local surface albedo by 31%–45%. Over the melting season, single ponding events increased the energy deposition of the sea ice by 35% compared to adjacent bare ice. Thus, single melt ponds may impact the summer energy budget as much as seasonal evolution over 1 month. Absorptivity and transmittance showed strong temporal and spatial variabilities independently of surface conditions, possibly due to the different internal sea-ice properties and under-ice biological processes. The differences in seasonal evolution shown for different sea-ice conditions strongly impacted the partitioning of shortwave solar radiation. This study shows that the formation and development of melt ponds, in reducing albedo by a third of bare ice sites, can notably increase the total summer heat deposition. The vastly different seasonal evolutions, different sea-ice conditions, and timing and duration of ponding events need to be considered when comparing local in-situ observations with large-scale satellite remote sensing datasets, which we suggest can help to improve numerical models.

Turbulent flow characteristics over rough permeable and impermeable gravel-bed stream- an experimental study

Abstract The study aims to presents flow heterogeneity over rough permeable and impermeable gravel-bed stream. The rough permeable stream is prepared by laying multiple layers of gravel, whereas its impermeable counterpart is presented by a resin-casted gravel-bed. In general, a common approach can be found in literature to mimic an impermeable bed by laying single layer of gravels, therefore some results are compared with single layer gravel-bed stream keeping the hydraulic conditions same. An acoustic Doppler velocimeter was used for flow measurements whereas; double averaging (DA) methodology was adopted for data analysis. The larger flow penetration depth and intense flow mixing in permeable gravel-bed infer sufficient impetus for organized flow turbulence and damping of DA Reynolds shear stresses whereas; the wall-blocking in resin-casted gravel-bed prevents fluid infiltration which leads to follow the linear stress profile away from the crest level. The damping of DA Reynolds shear stress (RSS) is compensated by enhanced DA form-induced shear stress (FISS). The results are further analysed under the light of the energy budget to characterize the mass-momentum exchange as it penetrates the subsurface layers. The energy budget indicates negative pressure energy diffusion rates corroborating gain in turbulence production in the permeable gravel-bed stream.

The Role of Turbulence Transport in Mechanical Energy Budgets

Abstract In this article, we study the role of turbulence transport on loss prediction using high fidelity scale-resolving simulations. For this purpose, we use eight high-fidelity simulation datasets and compute flux of entropy and stagnation pressure transport equation budgets for all the cases. We find that under certain unsteady inflow conditions, the stagnation pressure coefficient is not a reliable loss metric even at low Mach number conditions. This is due to the turbulence transport terms. The impact of these terms, typically assumed to be negligible, made stagnation pressure loss coefficient underpredict loss from 0% to over 40% when compared to entropy loss coefficient for cases considered here. This effect was most pronounced for the cases with highly unsteady inflow conditions and at low Reynolds numbers.

Debris cover effects on energy and mass balance of Batura Glacier in the Karakoram over the past 20 years

Abstract. The influence of supraglacial debris cover on glacier mass balance in the Karakoram is noteworthy. However, understanding of how debris cover affects the seasonal and long-term variations in glacier mass balance through alterations in the glacier's energy budget is incomplete. The present study coupled an energy–mass balance model with heat conduction within debris layers on debris-covered Batura Glacier in Hunza Valley to demonstrate the influence of debris cover on glacial surface energy and mass exchanges during 2000–2020. The mass balance of Batura Glacier is estimated to be -0.262±0.561 m w.e. yr−1, with debris cover accounting for a 45 % reduction in the negative mass balance. Due to the presence of debris cover, a significant portion of incoming energy is utilized for heating debris, leading to a large energy emission to the atmosphere via thermal radiation and turbulent sensible heat. This, in turn, reduces the melt latent heat energy at the glacier surface. We found that the mass balance exhibits a pronounced arch-shaped structure along the elevation gradient, which is associated with the distribution of debris thickness and the increasing impact of debris cover on the energy budget with decreasing elevation. Through a comprehensive analysis of the energy transfer within each debris layer, we have demonstrated that the primary impact of debris cover lies in its ability to modify the energy flux reaching the surface of the glacier. Thicker debris cover results in a smaller temperature gradient within debris layers, consequently reducing energy reaching the debris–ice interface. Over the past 2 decades, Batura Glacier has exhibited a trend towards less negative mass balance, likely linked to a decrease in air temperature and reduced ablation in areas with thin or sparse debris cover.

Towards user-friendly all-sky surface longwave downward radiation from space: General scheme and product

Abstract Longwave downward radiation (LWDR) is an important driving parameter in climate and hydrological models. Compared to traditional ground-based measurements, remote sensing has unique advantages in estimating global LWDR. However, for current remote sensing missions, as the typical available satellite-derived LWDR product with global coverage and hourly temporal resolution, the Clouds and the Earth’s Radiant Energy System-Synoptic (CERES-SYN) top of atmosphere and surface fluxes and clouds has a low spatial resolution (1° × 1°). There is still much room for improvement of the existing remote sensing LWDR products in terms of accuracy, spatiotemporal resolutions, and the ability to explain and quantify the changes of longwave radiation at various scales. To overcome these limitations, this paper developed a new global LWDR product with improved accuracy (the RMSE < 30 W/m2 over the globe), high temporal resolution (hourly), and spatial resolution (5 km) based on Moderate Resolution Imaging Spectroradiometer (MODIS) measurements. It serves as a LWDR product within the Long-term Earth System spatiotemporally Seamless Radiation budget dataset (referred to as LessRad). As the first long-term high-resolution, spatiotemporally continuous LWDR product (2002-2022, 1 h, 5 km), the LessRad reveals its advantages in studying the spatiotemporal variability of LWDR on finer scales. It also provides a valuable data source for various applications, such as analyzing land-atmosphere interactions and quantifying climate feedback, and thus, is potentially helpful for understanding the earth’s energy budget and dynamics.

Biophysical Drivers of Seasonal Hysteresis of Urban Heat Islands Across Climates and Urban Landscapes

AbstractThe phenomenon of seasonal hysteresis between urban heat islands (UHIs) and background surface temperature was reported in many studies, but a unifying explanation for the causes of its magnitude, shape and looping direction remained less clear, and the impacts of urban landscapes and background climate on such hysteresis were also the subject of enquiry. Here we performed a systematic study based on a long‐term offline simulation for surface energy budgets, surface air and skin temperature of the urban‐biosphere system in East China. It was carried out with an up‐to‐date land surface model — the CLM5‐LCZs (the Community Land Model version 5 (CLM5) coupled with a newly developed urban canopy model representing diverse urban landscapes with the local climate zones classification scheme). To isolate the causes of UHIs hysteresis, a simplified framework was constructed by combining a mathematical representation of grid/subgrid model outputs, the two‐resistance mechanism method and the time‐lagged cross correlation analysis. The results exhibited distinct hysteresis patterns of UHIs such as twisted, concave‐down, concave‐up and convex‐up curve across diverse background climates and urban landscapes. The basic shapes of hysteresis depict a gradual transition from concave/convex‐up to concave/convex‐down by characterized by a larger winter UHIs from a wetter and warmer climate to a drier and colder climate. Magnitude and seasonality of surface UHIs are separately dominated by the urban‐rural contrast of evapotranspiration and heat storage at daytime and nighttime, whereas across‐climate variations of daytime and nighttime surface UHIs hysteresis are regulated by the phase shift of convection efficiency and anthropogenic heat. Urban landscapes affect the hysteresis mainly by altering the amplitude rather than the shape of looping curves. Our work could provide a fuller picture of how to mitigate urban warming considering inter‐seasonal tradeoffs over a broader region.

Oscillatory behavior of drops in an ambient liquid under an impact

In this work, the oscillatory behavior of drops in an ambient liquid under an impact is numerically and experimentally studied using the Basilisk C and a self-designed falling drop tower facility, respectively. The numerical results of the oscillatory drops are consistent with the experimental data. A shape parameter e is proposed to describe the oscillatory behavior of the drops. Using the energy budgets of the drops, the influence of dimensionless parameters, i.e., Reynolds number Re, Ohnesorge number Oh, Bond number Bo, and viscosity ratio η on the oscillatory amplitude and period of the shape parameter e is systematically investigated. Furthermore, a scaling law between the maximum deformation of the drops during the oscillatory behavior and the dimensionless parameters is established, which can be interpreted by a mass-spring-damp model. It is found that the established model agrees well with our numerical and experimental results without any fitting parameters when the viscosity of the ambient liquid is near or smaller than that of the drop, and the maximum deformation has a clear proportionality as Oh2Re(1+Re).

Latitudinal variation in thermal performance of the common coral, Pocillopora spp.

Understanding how tropical corals respond to temperatures is important to evaluating their capacity to persist in a warmer future. We studied the common Pacific coral Pocillopora over 44° of latitude, and used populations at three islands with different thermal regimes to compare their responses to temperature using thermal performance curves (TPCs) for respiration and gross photosynthesis. Corals were sampled in the local autumn from Moorea, Guam, and Okinawa where mean (± s.d.) annual seawater temperature is 28.0±0.9°C, 28.9±0.7°C, and 25.1±3.4°C, respectively. TPCs for respiration were similar among latitudes, the thermal optimum (Topt) was above the local maximum temperature at all three islands, and maximum respiration was lowest at Okinawa. TPCs for gross photosynthesis were wider, implying greater thermal eurytopy, with a higher Topt in Moorea versus Guam and Okinawa. Topt was above the maximum temperature in Moorea, but was similar to daily temperatures over 13% of the year in Okinawa, and 53% of the year in Guam. There was greater annual variation in daily temperatures in Okinawa than Guam or Moorea, which translated to large variation in the supply of metabolic energy and photosynthetically fixed carbon at higher latitudes. Despite these trends, the differences in TPCs for Pocillopora were not profoundly different across latitudes, reducing the likelihood that populations of these corals could better match their phenotypes to future more extreme temperatures through migration. Any such response would place a premium on high metabolic plasticity and tolerance of large seasonal variations in energy budgets.

Similarity solutions in elastohydrodynamic bouncing

We investigate theoretically and numerically the impact of an elastic sphere on a rigid wall in a viscous fluid. Our focus is on the dynamics of the contact, employing the soft lubrication model in which the sphere is separated from the wall by a thin liquid film. In the limit of large sphere inertia, the sphere bounces and the dynamics is close to the Hertz theory. Remarkably, the film thickness separating the sphere from the wall exhibits non-trivial self-similar properties that vary during the spreading and retraction phases. Leveraging these self-similar properties, we establish the energy budget and predict the coefficient of restitution for the sphere. The general framework derived here opens many perspectives to study the lubrication film in impact problems.

Joint In-Orbit Computation and Communication for Minimizing Download Time From LEO Satellites

Downloading a large amount of data from a low Earth orbit satellite to a ground station can be challenging due to the limited contact window, dynamic channel quality, solar energy supply, and thermal management without an atmosphere. Considering such dynamics, this paper proposes a joint design of in-orbit computation and communication for download time minimization. We combine the non-convex thermal constraints and energy constraints into unified energy budget constraints with upper bound approximation, and computational efficiency is achieved by decomposing the resulting large-scale problem into a non-convex communication sub-problem, a convex computation sub-problem solvable with interior point method and a master problem that optimizes the energy budget allocation between computation and communication. The communication sub-problem is solved with a generalized-benders-decomposition-based algorithm that decouples downlink scheduling and power allocation based on a closed-form solution of optimal dual variables in the power allocation primal problem. And the master problem is solved with ternary search by proving the minimal download time is quasi-convex with respect to the energy budget allocation between computation and communication. Simulation results demonstrate that the proposed solution effectively reduces the download time, especially under strict energy constraints and severe channel variations.

Causes of Past African Temperature Change in PMIP Simulations of the Mid‐Holocene

AbstractCurrent‐generation climate models project that Africa will warm by up to 5°C in the coming century, severely stressing African populations. Past and ongoing work indicates, however, that the models used to create these projections do not match proxy records of past temperature in Africa during the mid‐Holocene (MH), raising concerns that their future projections may house large uncertainties. Rather than reproducing proxy‐based reconstructions of MH warming relative to the Pre‐Industrial (PI), models instead simulate MH temperatures very similar to or slightly colder than the PI. This data‐model mismatch could be due to a variety of factors, including biases in model surface energy budgets or inaccurate representation of the feedbacks between temperature and hydrologic change during the “Green Sahara.” We focus on the differences among model simulations in the Paleoclimate Modeling Intercomparison Project Phases 3 and 4 (PMIP3 and PMIP4), examining surface temperature and energy budgets to investigate controls on temperature and the potential model sources of this paleoclimate data‐model mismatch. Our results suggest that colder conditions simulated by PMIP3 and PMIP4 models during the MH are in large part due to the joint impacts of feedback uncertainties in response to increased precipitation, a strengthened West African Monsoon (WAM) in the Sahel, and the Green Sahara. We extend these insights into suggestions for model physics and boundary condition changes, and discuss implications for the accuracy of future climate model projections over Africa.

Pressure Balance and Energy Budget of the Nuclear Superbubble of NGC 3079

Superbubbles in the nuclear region of galaxies could be produced by active galactic nuclei (AGN) or nuclear starburst via different driving forces. We report analysis of the multiwavelength data of the kiloparsec-scale nuclear superbubble in NGC 3079 in order to probe the mechanisms driving the expansion of the superbubble. Based on the Chandra X-ray observations, we derive the hot-gas thermal pressure inside the bubble, which is about 1 order of magnitude higher than that of the warm ionized gas traced by optical lines. We derive a [C ii]-based star formation rate (SFR) of SFR ∼ 1.3 M ⊙ yr−1 from the nuclear region using the SOFIA/FIFI-LS observation. This SFR infers a radiation pressure toward the bubble shells much lower than the thermal pressure of the gases. The Very Large Array radio image infers that magnetic pressure at the northeast cap above the superbubble is less than the thermal pressure of the hot gas enclosed in the bubble but has a clearly larger extension. The magnetic field may thus still help to reconcile the expansion of the bubble. The observed thermal energy of the hot gas enclosed in the bubble requires an energy injection rate of ≳1042 erg s−1 within the bubble’s dynamical age, which is probably larger than the power provided by the current nuclear starburst and the parsec-scale jet. If this is true, stronger past AGN activity may provide an alternative energy source to drive the observed bubble expansion.

Pollution offsets the rapid evolution of increased heat tolerance in a natural population

Despite the increasing evidence for rapid thermal evolution in natural populations, evolutionary rescue under global warming may be constrained by the presence of other stressors. Highly relevant in our polluted planet, is the largely ignored evolutionary trade-off between heat tolerance and tolerance to pollutants. By using two subpopulations (separated 40 years in time) from a resurrected natural population of the water flea Daphnia magna that experienced a threefold increase in heat wave frequency during this period, we tested whether rapid evolution of heat tolerance resulted in reduced tolerance to the widespread metal zinc and whether this would affect heat tolerance upon exposure to the pollutant. Our results revealed rapid evolution of increased heat tolerance in the recent subpopulation. Notably, the sensitivity to the metal tended to be stronger (reduction in net energy budget) or was only present (reductions in heat tolerance and in sugar content) in the recent subpopulation. As a result, the rapidly evolved higher heat tolerance of the recent subpopulation was fully offset when exposed to zinc. Our results highlight that the many reports of evolutionary rescue to global change stressors may give a too optimistic view as our warming planet is polluted by metals and other pollutants.

Direct Numerical Simulation Database of High-Speed Flow over Parameterized Curved Walls

This study presents a direct numerical simulation (DNS) database of high-speed turbulent boundary layers (TBLs) subject to pressure gradients due to parametrically varied backward-facing and forward-facing wall curvatures, with an inflow Mach number of 4.9 and a friction Reynolds number of Reτ≃1000 immediately before the onset of wall curvature. The Mach and Reynolds numbers are significantly higher than those reported in the literature for the DNS of pressure-gradient TBLs. The flow conditions and baseline wall geometries are representative of experiments in the high-speed blowdown wind tunnel at the National Aerothermochemistry Laboratory at Texas A&M University. The wall steepness of the baseline geometry for both the backward-facing and forward-facing walls was systematically varied to cause attached, incipiently separated, and fully separated flows. Precomputed flow statistics, including turbulent kinetic energy budgets, are available on the website of the Turbulence Modeling Resource of the NASA Langley Research Center, allowing other investigators to query any property of interest.

Axial friction coefficient of turbulent spiral Poiseuille flows

Direct numerical simulations of spiral Poiseuille flows in a narrow gap geometry are performed with the aim of identifying the mechanisms governing the dynamics of the axial friction coefficient. The investigation has explored a small portion of the Reynolds number–Taylor number phase space ( $600 \leq Re \leq 5766$ and $1500 \leq Ta \leq 5000$ ), for which reference experimental results are available. The study is focused on the mechanism leading to the enhancement of the axial friction coefficient with the Taylor number when the Reynolds number is kept constant. The analysis of the spatial distribution of the Reynolds stress tensor and of the turbulent energy budget has evidenced the key role of the pressure–strain correlation in the energy transfer from the azimuthal to the axial component. The latter eventually determines the increase of the axial friction coefficient through the enhanced radial mixing of axial momentum. Data have also shown that the flow dynamics is heavily dependent on the $Ta/Re$ ratio, and different regimes develop (ranging from laminar to turbulent), each with peculiar behaviours.

From snow accumulation to snow depth distributions by quantifying meteoric ice fractions in the Weddell Sea

Abstract. Year-round snow cover is a characteristic of the entire Antarctic sea ice cover, which has significant implications for the energy and mass budgets of sea ice, e.g., by preventing surface melt in summer and enhancing sea ice growth through extensive snow ice formation. However, substantial observational gaps in the seasonal cycle of Antarctic sea ice and its snow cover limit the understanding of important processes in the ice-covered Southern Ocean. They also introduce large uncertainties in satellite remote sensing applications and climate studies. Here we present results from 10 years of autonomous snow observations from Snow Buoys in the Weddell Sea. To distinguish between actual snow depth and potential snow ice thickness within the accumulated snowpack, a one-dimensional thermodynamic sea ice model is applied along the drift trajectories of the buoys. The results show that potential snow ice formation, with an average maximum thickness of 35 cm, was detected along 41 % of the total track length of the analyzed Snow Buoy tracks, which corresponds to about one-quarter of the snow accumulation. In addition, we simulate the evolution of internal snow properties along the drift trajectories with the more complex SNOWPACK model, which results in superimposed ice thicknesses between 0 and 14 cm on top of the snow ice layer. These estimates will provide an important reference dataset for both snow depth and meteoric ice rates in the Southern Ocean.

Energy Budget Characterisation of the Optimal Disturbance in Stratified Shear Flow

Stratified Taylor–Couette flow (STCF) undergoes transient growth. Recent studies have shown that there exists transient amplification in the linear regime of counter-rotating STCF. The kinetic budget of the optimal transient perturbation is analysed numerically to simulate the interaction of the shear production (SP), buoyancy flux (BP), and other energy components that contributes to the total optimal transient kinetic energy. These contributions affect the total energy by influencing the perturbation to extract kinetic energy (KE) from the mean flow. The decay of the amplification factor resulted from the positive amplification of both BP and SP, while the growth is attributed to the negative and positive amplification of BP and SP, respectively. The optimal SP is positively amplified, implying that there is the possibility of constant linear growth. These findings agree with the linear growth rate for increasing values of Grashof number.

Dynamics of Carbon and Water Fluxes over Cropland and Agroforest Ecosystems on the Southern Chinese Loess Plateau

Studies on the spatiotemporal dynamics in ecosystem carbon and water exchanges are essential in predicting the effects of climate change on regional carbon and energy budgets. Using the eddy covariance technique, carbon and water fluxes were observed in a typical winter wheat ecosystem (WWE) and an agroforest ecosystem (AFE) in the southern Loess Plateau from 2004 to 2010. The seasonal and inter-annual variability in gross primary productivity (GPP), net ecosystem exchange (NEE), evapotranspiration (ET), and water use efficiency (WUE) were examined and the main influencing factors were identified using the Pearson correlation. The results indicate that the seasonal GPP and NEE showed a bimodal distribution in WWE, while this was unimodal in AFE. The sinusoidal function did well in the characterization of seasonal ET dynamics for both ecosystems, with the determination coefficients being 0.85 and 0.94, respectively. In WWE and AFE, the annual mean GPP were 724.33 and 723.08 g C m−2 a−1, respectively, and the corresponding ET were 392.22 and 410.02 mm a−1. However, the difference in NEE between the two ecosystems was obvious, NEE were −446.28 and −549.08 g C m−2 a−1, respectively, showing a stronger carbon sink in AFE. There were strong coupling relationships between the GPP and ET of both ecosystems; the overall slopes were 1.71 and 1.69, respectively. The seasonal trend of WUE was bimodal in WWE, with peak values of 3.94 and 3.65 g C kg−1 H2O, occurring in November and April, respectively. However, the monthly WUE in AFE had one single peak of 4.07 g C kg−1 H2O in January. Photosynthetically active radiation (PAR) and soil temperature (Ts) were most positively correlated with GPP, net radiation (Rn) and Ts were the major factors influencing ET, while vapor pressure deficit (VPD) and soil water content (SWC) were the major influencing factors for WUE. These results provide observational support for regional carbon neutrality simulations.

Finite Reynolds Number Effect on Small-Scale Statistics of Homogeneous Isotropic Turbulence

Since about 1997, the realisation that the finite Reynolds number (FRN) effect needs to be carefully taken into account when assessing the behaviour of small-scale statistics came to the fore. The FRN effect can be analysed either in the real domain or in the spectral domain via the scale-by-scale energy budget equation or the transport equation for the energy spectrum. This analysis indicates that the inertial range (IR) is established only when the Taylor microscale Reynolds number Reλ is infinitely large, thus raising doubts about published power-law exponents at finite values of Reλ, for either the second-order velocity structure function (δu)2¯ or the energy spectrum. Here, we focus on the transport equation of (δu)2¯ in decaying grid turbulence, which represents a close approximation to homogeneous isotropic turbulence. Regarding small-scale effects, the large-scale forcing term associated with streamwise advection decreases as Reλ increases and finally disappears when Reλ is sufficiently large. An approach based on the dual scaling of (δu)2¯, i.e., a scaling based on the Kolmogorov scales (when the separation r is small) and another based on the integral scales (when r is large), yields (δu)2¯∼r2/3 when Reλ is infinitely large. This approach also yields (δu)n¯∼rn/3 when Reλ is infinitely large. These results seem to be supported by the trend as Reλ increases according to the available experimental data. Overall, the results for decaying turbulence strongly suggest that a tendency towards the predictions of K41 cannot be dismissed at least at Reynolds numbers that are currently beyond the reach of experiments and direct numerical simulations.

The ARGOS Instrument for Stratospheric Aerosol Measurements

Atmospheric aerosols represent an important component of the Earth’s climate system because they can contribute both positive and negative forcing to the energy budget. We are developing the Aerosol Radiometer for Global Observations of the Stratosphere (ARGOS) instrument to provide improved measurements of stratospheric aerosols in a compact package. ARGOS makes limb scattering measurements from space in eight directions simultaneously, using two near-IR wavelengths for each viewing direction. The combination of forward and backward scattering views along the orbit track gives additional information to constrain the aerosol phase function and size distribution. Cross-track views provide expanded spatial coverage. ARGOS will have a demonstration flight through a hosted payload provider in the fall of 2024. The instrument has completed pre-launch environmental testing and radiometric characterization tests. The hosted payload approach offers advantages in size, weight, and power margins for instrument design compared to other approaches, with significant benefits in terms of reducing infrastructure requirements for the instrument team.